* IV REGRESSION OF EFFECT OF EXPECTED INFLATION 12 MONTHS AHEAD ON:
*		OWN PRICE ANNUAL CHANGE - EXPECTED AT TIME OF INTERVIEW, OBSERVED AT T+1, T+2, T+4
*		LOG LABOR DEMAND  - OBSERVED AT T+1, T+2, T+3, T+4


* ESTIMATE SAME MODEL IN 4 WAYS: WITH AND WITHOUT TRIMMING DATA AND WITH AND WITHOUT SAMPLE WEIGHTS.
* POINT RESULTS OF ALL COMBINATIONS (WEIGHT/TRIM) ARE REPORTED IN TABLE G1-G8.
* FIG 7 PLOTS THE CDF OF TSTATS OF ESTIMATED EFFECTS ON OWN PRICE CHANGES AND LABOR DEMAND FROM UNWEIGHTED REGRESSIONS ON TRIMMED SAMPLE.

* FIRST LOOP (OVER WEIGHTED/UNWEIGHTED) ESTIMATES TRIMMED WEIGHTED AND TRIMMED UNWEIGHTED RESULTS
* SECOND LOOP (OVER WEIGHTED/UNWEIGHTED) ESTIMATES UNTRIMMED WEIGHTED AND UNTRIMMED UNWEIGHTED RESULTS

* COEFFICIENTS AND SES ARE STORED IN MATRIXES (DP - OWN PRICE; LD - LAB DEMAND)
* MATRIXES ARE THEN SAVED AS DTA AND USED TO GENERATE TABLES AND FIGURES

* TRIMMED WEIGHTED/UNWEIGHTED IV

quietly{

use WorkingData, clear

forvalues w=0/1{
	
cap drop mypeso
if $_w==1{
gen mypeso = peso /* weighted est */
}
if $_w==0{
gen mypeso = 1    /* unweighted est */
}


local explvar=" dpre i.dim i.area i.setton5 "


* OWN PRICES, expected, obs next quarter, intwo quarter in 4 quarters

mat def DP=J($ntrims, 9, .)
local r=1
forvalues q=$firsttrim/$lasttrim{
mat DP[$_r,1 ] = $_q

local c=2
local cc=3
foreach y in dprez F1DPZ F2DPZ F4DPZ{

cap qui ivregress 2sls $_y $_explvar (it12all=treatment) [aw=mypeso] if trims !=1 & dp==$_q, vce(robust)
if _rc==0{
mat DP[$_r,$_c ] = _b[it12all]
mat DP[$_r,$_cc] = _se[it12all]
}

local c  = $_c +2
local cc = $_cc+2

}

local r = $_r + 1
}


 * (log) Labor demand 1, 2, 3, 4 quarers after info assignment


mat def LD=J($ntrims, 9, .)
local r=1
forvalues q=$firsttrim/$lasttrim{
mat LD[$_r,1 ] = $_q

local c=2
local cc=3
foreach y in F1LAD F2LAD F3LAD F4LAD{

cap qui ivregress 2sls $_y $_explvar (it12all=treatment) [aw=mypeso] if trims !=1 & dp==$_q, vce(robust)
if _rc==0{
mat LD[$_r,$_c ] = _b[it12all]
mat LD[$_r,$_cc] = _se[it12all]
}

local c  = $_c +2
local cc = $_cc+2

}

local r = $_r + 1
}


* Store coefficient estiamtes and robust SE in dtas 

* own prices 
preserve
drop _all
svmat DP
rename DP1 dp
format dp %tq

local i=2
local ii=3
foreach y in dprez F1DPZ F2DPZ F4DPZ{
rename DP$_i bb_$_y
rename DP$_ii se_$_y

local i = $_i + 2
local ii = $_ii + 2
}


save DataTabIV2SLS_DPre_W$_w, replace

restore

* labor demand

preserve
drop _all
svmat LD
rename LD1 dp
format dp %tq

local i=2
local ii=3
foreach y in F1LAD F2LAD F3LAD F4LAD {
rename LD$_i bb_$_y
rename LD$_ii se_$_y

local i = $_i + 2
local ii = $_ii + 2
}


save DataTabIV2SLS_LabDem_W$_w, replace

restore

}

}

* UNTRIMMED WEIGHTED/UNWEIGHTED IV

quietly{

use WorkingData, clear

forvalues w=0/1{
	
cap drop mypeso
if $_w==1{
gen mypeso = peso /* weighted est */
}
if $_w==0{
gen mypeso = 1    /* unweighted est */
}


local explvar=" dpre i.dim i.area i.setton5 "


* OWN PRICES, expected, obs next quarter, intwo quarter in 4 quarters

mat def DP=J($ntrims, 9, .)
local r=1
forvalues q=$firsttrim/$lasttrim{
mat DP[$_r,1 ] = $_q

local c=2
local cc=3
foreach y in dprez F1DPZ F2DPZ F4DPZ{

cap qui ivregress 2sls $_y $_explvar (it12all=treatment) [aw=mypeso] if dp==$_q, vce(robust)
if _rc==0{
mat DP[$_r,$_c ] = _b[it12all]
mat DP[$_r,$_cc] = _se[it12all]
}

local c  = $_c +2
local cc = $_cc+2

}

local r = $_r + 1
}


* (log) labor demand 1, 2, 3, 4 qtrs after info assignemnt


mat def LD=J($ntrims, 9, .)
local r=1
forvalues q=$firsttrim/$lasttrim{
mat LD[$_r,1 ] = $_q

local c=2
local cc=3
foreach y in F1LAD F2LAD F3LAD F4LAD{

cap qui ivregress 2sls $_y $_explvar (it12all=treatment) [aw=mypeso] if  dp==$_q, vce(robust)
if _rc==0{
mat LD[$_r,$_c ] = _b[it12all]
mat LD[$_r,$_cc] = _se[it12all]
}

local c  = $_c +2
local cc = $_cc+2

}

local r = $_r + 1
}


* store coeff estimates and robust SE in dtas

* own prices 
preserve
drop _all
svmat DP
rename DP1 dp
format dp %tq

local i=2
local ii=3
foreach y in dprez F1DPZ F2DPZ F4DPZ{
rename DP$_i bb_$_y
rename DP$_ii se_$_y

local i = $_i + 2
local ii = $_ii + 2
}


save DataTabIV2SLS_U_DPre_W$_w, replace

restore

* labor demand

preserve
drop _all
svmat LD
rename LD1 dp
format dp %tq

local i=2
local ii=3
foreach y in F1LAD F2LAD F3LAD F4LAD {
rename LD$_i bb_$_y
rename LD$_ii se_$_y

local i = $_i + 2
local ii = $_ii + 2
}


save DataTabIV2SLS_U_LabDem_W$_w, replace

restore

}

	
	
}

* *************************** GENERATE TABLES G1-G8 *************************** 

* ******************************** OWN PRICES *********************************

* TRIMMED UNWEIGHTED
use DataTabIV2SLS_DPre_W0, clear

quietly{
foreach y in dprez F1DPZ F2DPZ F4DPZ{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_Dpre_W0, replace
}

format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll
sort dp
outsheet dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll using TabIV2SLSCoef_Dpre_W0, replace delim("&")
}

* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* DPREZ - EXPECTED OWN PRICE CHANGE IN 12 MONTHS
* FXDPZ - OBSERVED ANNUAL PRICE CHANGE XQTR AFTER ASSINGMENT
sort dp 
* *******************************************
* *** TAB G1: TRIMMED UNWEIGHTED IV REGS ***
* *******************************************
list dp bb_dprez sig_dprez se_dprez bb_F1DPZ sig_F1DPZ se_F1DPZ bb_F2DPZ sig_F2DPZ se_F2DPZ bb_F4DPZ sig_F4DPZ se_F4DPZ , linesize(255) noobs header
* =============================================================================


* UNTRIMMED UNWEIGHTED 
use DataTabIV2SLS_U_DPre_W0, clear

quietly{
foreach y in dprez F1DPZ F2DPZ F4DPZ{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_U_Dpre_W0, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll
sort dp
outsheet dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll using TabIV2SLSCoef_U_Dpre_W0, replace delim("&")
}

* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* DPREZ - EXPECTED OWN PRICE CHANGE IN 12 MONTHS
* FXDPZ - OBSERVED ANNUAL PRICE CHANGE XQTR AFTER ASSINGMENT
sort dp
* ************************************************
* **** TAB G3: UNTRIMMED UNWEIGHTED IV REGS ****
* ************************************************
list dp bb_dprez sig_dprez se_dprez bb_F1DPZ sig_F1DPZ se_F1DPZ bb_F2DPZ sig_F2DPZ se_F2DPZ bb_F4DPZ sig_F4DPZ se_F4DPZ , linesize(255) header noobs
* =============================================================================


* TRIMMED WEIGHTED 
use DataTabIV2SLS_DPre_W1, clear

quietly{
foreach y in dprez F1DPZ F2DPZ F4DPZ{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_Dpre_W1, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll
sort dp
outsheet dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll using TabIV2SLSCoef_Dpre_W1, replace delim("&")
}

* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* DPREZ - EXPECTED OWN PRICE CHANGE IN 12 MONTHS
* FXDPZ - OBSERVED ANNUAL PRICE CHANGE XQTR AFTER ASSINGMENT
sort dp
* ******************************************
* **** TAB G5: TRIMMED WEIGHTED IV REGS ****
* ******************************************
list dp bb_dprez sig_dprez se_dprez bb_F1DPZ sig_F1DPZ se_F1DPZ bb_F2DPZ sig_F2DPZ se_F2DPZ bb_F4DPZ sig_F4DPZ se_F4DPZ , linesize(255) header noobs
* =============================================================================


* UNTRIMMED WEIGHTED 
use DataTabIV2SLS_U_DPre_W1, clear

quietly{
foreach y in dprez F1DPZ F2DPZ F4DPZ{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_U_Dpre_W1, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll
sort dp
outsheet dp bb_dprez sig_dprez se_dprez z1 bb_F1DPZ sig_F1DPZ se_F1DPZ z2 bb_F2DPZ sig_F2DPZ se_F2DPZ z3 bb_F4DPZ sig_F4DPZ se_F4DPZ ll using TabIV2SLSCoef_U_Dpre_W1, replace delim("&")

}

* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* DPREZ - EXPECTED OWN PRICE CHANGE IN 12 MONTHS
* FXDPZ - OBSERVED ANNUAL PRICE CHANGE XQTR AFTER ASSINGMENT
sort dp 
* **********************************************
* * **** TAB G7: UNTRIMMED WEIGHTED IV REGS ****
* **********************************************
list dp bb_dprez sig_dprez se_dprez bb_F1DPZ sig_F1DPZ se_F1DPZ bb_F2DPZ sig_F2DPZ se_F2DPZ bb_F4DPZ sig_F4DPZ se_F4DPZ  , linesize(255) header noobs
* =============================================================================


* ********************* LABOR DEMAND *********************************************************


* TRIMMED UNWEIGHTED REG 
use DataTabIV2SLS_LabDem_W0, clear

quietly{
foreach y in F1LAD F2LAD F3LAD F4LAD{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_LabDem_W0, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq


* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll
sort dp
outsheet dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll using TabIV2SLSCoef_LabDem_W0, replace delim("&")
}


* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* FXLAD - LABOR DEMAND (LOG) X QTR AFTER ASSINGMENT
sort dp
* *********************************************
* **** TAB G2: TRIMMED UNWEIGHTED IV REGS ****
* *********************************************
list dp bb_F1LAD sig_F1LAD se_F1LAD bb_F2LAD sig_F2LAD se_F2LAD bb_F3LAD sig_F3LAD se_F3LAD bb_F4LAD sig_F4LAD se_F4LAD  , linesize(255) header noobs
* =============================================================================


* UNTRIMMED UNWEIGHTED REG
use DataTabIV2SLS_U_LabDem_W0, clear

quietly{
foreach y in F1LAD F2LAD F3LAD F4LAD{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_U_LabDem_W0, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll
sort dp
outsheet dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll using TabIV2SLSCoef_U_LabDem_W0, replace delim("&")
}


* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* FXLAD - LABOR DEMAND (LOG) X QTR AFTER ASSINGMENT
sort dp 
* **********************************************
* **** TAB G4: UNTRIMMED UNWEIGHTED IV REGS ****
* **********************************************
list dp bb_F1LAD sig_F1LAD se_F1LAD bb_F2LAD sig_F2LAD se_F2LAD bb_F3LAD sig_F3LAD se_F3LAD bb_F4LAD sig_F4LAD se_F4LAD  , linesize(255) header noobs
* =============================================================================


* TRIMMED WEIGHTED REG 
use DataTabIV2SLS_LabDem_W1, clear

quietly{
foreach y in F1LAD F2LAD F3LAD F4LAD{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_LabDem_W1, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll
sort dp
outsheet dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll using TabIV2SLSCoef_LabDem_W1, replace delim("&")
}

* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* FXLAD - LABOR DEMAND (LOG) X QTR AFTER ASSINGMENT
sort dp 
* ******************************************
* **** TAB G6: TRIMMED WEIGHTED IV REGS ****
* ******************************************
list dp bb_F1LAD sig_F1LAD se_F1LAD bb_F2LAD sig_F2LAD se_F2LAD bb_F3LAD sig_F3LAD se_F3LAD bb_F4LAD sig_F4LAD se_F4LAD  , linesize(255) header noobs
* =============================================================================


* UNTRIMMED WEIGHTED REG 
use DataTabIV2SLS_U_LabDem_W1, clear

quietly{
foreach y in F1LAD F2LAD F3LAD F4LAD{

gen p_$_y = 2*(1-normal(abs(bb_$_y / se_$_y)))
gen sig_$_y=""
replace sig_$_y="$^{*}$" if p_$_y<0.1
replace sig_$_y="$^{**}$" if p_$_y<0.05 
replace sig_$_y="$^{***}$" if p_$_y<0.01
save TableIV2SLS_U_LabDem_W1, replace
}



format bb_* %3.2f
format se_* %6.3f
format dp %tq

* these are inserted to automate subsequent editing and formatting of the table
gen z1="99"
gen z2="99"
gen z3="99"
gen ll="\\"
order dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll
sort dp
outsheet dp bb_F1LAD sig_F1LAD se_F1LAD z1 bb_F2LAD sig_F2LAD se_F2LAD z2 bb_F3LAD sig_F3LAD se_F3LAD z3 bb_F4LAD sig_F4LAD se_F4LAD ll using TabIV2SLSCoef_U_LabDem_W1, replace delim("&")
}

* =============================================================================
* PRINT TO LOG COEFFS (bb_), STATSIG (sig_) AND SES (se_)
* FXLAD - LABOR DEMAND (LOG) X QTR AFTER ASSINGMENT
sort dp 
* ********************************************
* **** TAB G8: UNTRIMMED WEIGHTED IV REGS ****
* ********************************************
list dp bb_F1LAD sig_F1LAD se_F1LAD bb_F2LAD sig_F2LAD se_F2LAD bb_F3LAD sig_F3LAD se_F3LAD bb_F4LAD sig_F4LAD se_F4LAD  , linesize(255) header noobs
* =============================================================================


* *************************** GENERATE TABLES G1-G8: END *************************** 

* *************************** GENERATE FIG 7 *************************** 

use DataTabIV2SLS_DPre_W0, clear

quietly{
foreach x in dprez F1DPZ F2DPZ F4DPZ{
gen tstat_$_x= bb_$_x / se_$_x
}
reshape long bb_ se_ tstat_, i(dp) j(v, string)
drop if tst==.

gen x=1
sort tsta
gen c=sum(x)
keep c tsta
gen v=1
save ForFigTstatIV2SLS, replace

use DataTabIV2SLS_LabDem_W0, clear

foreach x in F1LAD F2LAD F3LAD F4LAD{
gen tstat_$_x= bb_$_x / se_$_x
}
reshape long bb_ se_ tstat_, i(dp) j(v, string)
drop if tst==.

gen x=1
sort tsta
gen c=sum(x)
format tstat %3.1f
keep tstat c
gen v=2
append using ForFigTstatIV2SLS
save ForFigTstatIV2SLS, replace


lab def lv 1 "Expected and observed price changes" 2 "Observed labour demand"
lab val v lv
scatter c tsta, msy(oh) xli(-1.96 +1.96) lw(medthick) lc(red) yla(0(10)160, angle(0)) xla(-2(0.1)2, angle(90) labsi(*0.8)) xti("Tstat of IV estimate") yti("Count") xsize(10) ysize(4) by(v, note("") row(1) graphregion(c(white)) ) subti(,bc(white))

*graph export PostAER/FigTstatPricesLD_AER.eps, replace
graph export Fig7.eps, replace
}

* =============================================================================
* FIG 7 CDF of TSTATS OF CONVENTIONAL IV ESTIMATES FOR OWN PRICES AND LABOR DEMAND 
* **********************************
* ******* DATA FOR FIG 7 ***********
* **********************************
sort v c 
list v c tstat , linesize(255) header noobs
* =============================================================================



* *************************** GENERATE FIG 7: END *************************** 
